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Abstract 

We calculate the cross section of J/ij: inclusive production in lepton-nucleon 
deep-inelastic scattering via the weak charged current within the factorization for- 
malism of nonrelativistic quantum chromodynamics. The hadronic system pro- 
duced in association with the J/ij: meson can be of light-quark origin or contain 
charm. We take into account both direct production and feed-down from directly- 
produced heavier charmonia. We provide the cross sections of all contributing par- 
tonic subprocesses in analytic form. We present theoretical predictions for the J/il) 
transverse-momentum and rapidity distributions, which can be measured in current 
neutrino-nucleon scattering experiments at CERN and Fermilab, and possibly at 
the future electron-proton collider THERA at DESY. We find the cross section to 
be at the edge of observability at DESY HERA. 

PACS numbers: 12.38.-t, 13.60.Hb, 13.60.Le, 14.40.Gx 



1 Introduction 

Since the discovery of the J /ip meson in 1974, charmonium has provided a useful labora- 
tory for quantitative tests of quantum chromodynamics (QCD) and, in particular, of the 
interplay of perturbative and nonperturbative phenomena. The factorization formalism 
of nonrelativistic QCD (NRQCD) [1] provides a rigorous theoretical framework for the 
description of heavy- quarkonium production and decay. This formalism implies a separa- 
tion of short-distance coefficients, which can be calculated perturbatively as expansions in 
the strong-coupling constant a^, from long-distance matrix elements (MEs), which must 
be extracted from experiment. The relative importance of the latter can be estimated by 
means of velocity scaling rules, i.e., the MEs are predicted to scale with a definite power of 
the heavy-quark {Q) velocity v in the limit v <^1. In this way, the theoretical predictions 
arc organized as double expansions in Us and v. A crucial feature of this formalism is that 
it takes into account the complete structure of the QQ Fock space, which is spanned by 
the states n — ^s+i^W ^j^jj definite spin 5", orbital angular momentum L, total angular 
momentum J, and colour multiplicity C = 1, 8. In particular, this formalism predicts 
the existence of colour-octet (CO) processes in nature. This means that QQ pairs are 
produced at short distances in CO states and subsequently evolve into physical, colour- 
singlet (CS) quarkonia by the nonperturbative emission of soft gluons. In the limit w — > 0, 
the traditional CS model (CSM) [2] is recovered. The greatest triumph of this formahsm 
was that it was able to correctly describe [3] the cross section of inclusive charmonium 
hadroproduction measured in pp collisions at the Fermilab Tevatron [4] , which had turned 
out to be more than one order of magnitude in excess of the theoretical prediction based 
on the CSM. 

Apart from this phenomenological drawback, the CSM also suffers from severe concep- 
tual problems indicating that it is incomplete. These include the presence of logarithmic 
infrared divergences in the 0{as) corrections to P-wave decays to light hadrons and in 
the relativistic corrections to S'-wave annihilation [5], and the lack of a general argument 
for its validity in higher orders of perturbation theory. While the /c-r-factorization [6] and 
hard-comover-scattering [7] approaches manage to bring the CSM prediction much closer 
to the Tevatron data, they do not cure the conceptual defects of the CSM. The colour 
evaporation model [8], which is intuitive and useful for qualitative studies, also leads to 
a significantly better description of the Tevatron data, but it is not meant to represent 
a rigorous framework for perturbation theory. In this sense, a coequal alternative to the 
NRQCD factorization formalism is presently not available. 

In order to convincingly establish the phenomenological significance of the CO pro- 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction, ep and vN deep-inelastic 
scattering (DIS), e+e^ annihilation, 77 collisions, and 6-hadron decays may be found in 
the hterature; see Ref. [9] and references cited therein. Furthermore, the polarization of 
■0' mesons produced directly and of J/%1) mesons produced promptly, i.e., either directly 
or via the feed-down from heavier charmonia, which also provides a sensitive probe of 
CO processes, was investigated [10-13]. Until very recently, none of these studies was 
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able to prove or disprove the NRQCD factorization hypothesis. However, HI data of 

e + p^ e + J/ip + X in DIS at HERA [14] and DELPHI data of 7 + 7 ^ J/ip + X at 
LEP2 [15] provide first independent evidence for it [16,17]. 

Recently, we studied J/ip inclusive production in ep DIS via the electromagnetic cur- 
rent [16] and in uN DIS via the weak neutral current (NC) [18]. In both cases, we 
allowed for an additional hadron jet j in the final state, i.e., we considered the processes 
e + p—>-e + J/ip+j + X and u + N ^ p + J/ib + j + X, respectively, where X denotes 
the nucleus remnant. In this way, the inelasticity variable z, which measures the fraction 
of energy transferred from the virtual photon (7*) or Z boson (Z*) to the J/ijj meson 
in the nucleus rest frame, can take values below unity, away from the endpoint z — 1, 
where diffractive production takes place. At the same time, the J/ip meson can acquire 
finite transverse momentum p^ in the 7*p and Z*N centre-of-mass (CM) frames, and the 
hadronic system X' consisting of j and X can acquire finite mass Mx>. By the same to- 
ken, diffractive events can be eliminated from the experimental data sample by applying 
appropriate acceptance cuts on z, p^, or Mx'- 

In this paper, we complete our study of J/%1) inclusive production in DIS by considering 
the processes e + p v + J/ip + j + X and u -\- N — > e + J/il) + j + X , which are 
mediated by the weak charged current (CC). The former process can, in principle, be 
analyzed experimentally at the hadron-electron ring accelerator (HERA), with CM energy 
= 318 GeV, in operation at DESY, and at THERA [19], a future hadron-electron 
supercollider with y/S — 1 TeV that uses the HERA proton beam and the electron 
beam of the TcV-Energy Superconducting Linear Accelerator (TESLA), which is presently 
being designed and planned at DESY. The latter process can be investigated by the 
NuTeV Collaboration at Femilab (Experiment E815) [20] and by the CHORUS [21] and 
NOMAD [22] Collaborations at CERN. Since these experiments utilize i/^ or 17^ beams and 
J/ip mesons are preferably identified through their decays to e"^e~ or pairs, signal 

events would contain spectacular fi^e~^e~ or ^^fi~^fi~ trileptons together with hadronic 
debris from the target. In fact, events of both the former [23] and latter [24] types were 
observed in i/N fixed-target experiments at CERN and Fermilab. A theoretical study of 
the processes + N — > n~ + J/ijj + X and V^j, + N /j,'^ + J/i/j + X may be found 
in Ref. [25], where the J/ip mesons are assumed to arise electronically, i.e., from virtual 
photons that are emitted from charged leptons, quarks, or W bosons involved in the hard 
scatterings. To our knowledge, hadronic production mechanisms in the framework of 
NRQCD or the CSM have not yet been considered in the literature. On the other hand, 
open CC production in uN CC DIS was studied in QCD [26] . 

Diffractive contributions to prompt J/t/j production in CC DIS should be greatly 
suppressed. In contrast to the photon and Z boson, the W boson, being charged, cannot 
directly fluctuate into charmonium. A secondary CC interaction is necessary for that. 
We refer the interested reader to the pioneering literature [27]. An alternative set of 
Feynman diagrams that allows charmonium to diffractively scatter off the hadron target 
involves a timelike W* boson, which decays. The charmonium state would then originate 
from a 7* or Z* boson that is coupled to the lepton line, the W* boson, or the decay 
products of the latter. Both sets of Feynman diagrams interfere, and the resulting cross 
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section is suppressed by two powers of Sommerfeld's fine-structure constant a relative to 
its counterparts in electromagnetic or NC DIS. For this reason, diffractive processes are 
disregarded in the following. In particular, there is no need for acceptance cuts on z, p^, 
or Mx'. 

The J/ip mesons can be produced directly; or via radiative or hadronic decays of 
heavier charmonia, such as XcJ ^-nd ip' mesons; or via weak decays of b hadrons. The 
respective decay branching fractions are B{xco ^ J/ip + 7) = (1.02 ± 0.17)%, B{xci 
J/^ + 7) = (31.6 ± 3.2)%, B{xc2 ^ W + 7) = (18.7 ± 2.0)%, 5(^' ^ J/^ + X) = 
(55.7±2.6)%, and B{B J/ip + X) = (1.15±0.06)% [28]. The b hadrons can be detected 
by looking for displaced decay vertices with dedicated vertex detectors. Therefore and 
because of the smallness of B{B — > J/t/j + X), J/t/j production through 6-hadron decay is 
not considered here. The cross sections of the four residual indirect production channels 
may be approximated by multiplying the direct-production cross sections of the respective 
intermediate charmonia with their decay branching fractions to J/ijj mesons. 

This paper is organized as follows. In Section 2, we list our formulas for the cross 
sections of prompt J/t/j production in CC DIS. Lengthy expressions are relegated to the 
Appendix. In Section 3, we present our numerical results and discuss their phenomeno- 
logical implications. Our conclusions are summarized in Section 4. 



2 Analytic results 

In this section, we present our analytic results for the cross sections of 



Pi + N 

r + iv 



■ r + H+j + x, 

■l+ + H + j + X, 
■ui + H + j+X, 
■Vi + H + j + X, 



(1) 
(2) 
(3) 
(4) 



in CC DIS, where I — e, ji^T, N is the hadron target, and H — J/ip, XcJi i^' with J = 0, 1, 2. 
We work in the fixed-fiavour- number scheme, i.e., we have n/ = 3 active quark fiavours 

g = -u, d, s in the nucleus A^. We treat the nucleus N and the quarks q as massless. The 
charm quark c and antiquark c, with mass mc, only appear in the final state. The cc Fock 
states contributing at lowest order (LO) in v are n = ^S[^\ ^Sq^\ ^S[^\ ^Pj'^ ii H = J/ijj, ip' 
and n = '^Pf\'^Sf^ ii H = Xcj- Their MEs satisfy the multiplicity relations 



QXcJ 
QXcJ 



3p(8) 



3p(l) 
"35.(8) 



(2J + 1)(C»^ 

(2J + 1)(C^=« 
(2 J + 1) 



3p(8) 

3o(l) 



3^(8) 



(5) 



which follow to LO in v from heavy-quark spin symmetry. 
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We are thus led to consider the following partonic subprocesses 

z/; + a — > /~ + cc[n] + b, (6) 

Vi + a^t + cc[n] + b, (7) 

r + a ^ i/z + cc[n] + b, (8) 

+ a ^ 77; + cc[n] + b, (9) 

where (a, b) = {d, u), {u, d), {s. c) for the transitions diagonal w.r.t. the Cabibbo-Kobayashi- 
Maskawa (CKM) flavour mixing and (a, 6) = {s,u), {u,s), {d,c) for the CKM-suppressed 
ones. Here, it is understood that q — q. Typical Feynman diagrams are depicted in 
Figs. 1 (a) and (b). If b = c [see Fig. 1 (b)], then the jet j contains a charmed hadron, 
which can be identified experimentally. In this case, all the cc Fock states enumerated 
above contribute. In particular, this production mode is also possible within the CSM. 
By contrast, ifb = u,d,s [see Fig. 1 (a)], then only the CO channel n = ^S[^^ contributes, 
i.e., the prompt production of a J/ip meson in association with an otherwise charmless 
hadronic system in CC DIS is a pure CO process. 

We now define the kinematics, considering process (1) for definiteness. As in Refs. [16,18 
we denote the four-momenta of ui, N, l~ , H, and j by k, P, k', pn, and p', respectively. 
The parton struck by the virtual W boson (W*) carries four-momentum p = xP. We 
take the mass M of to be M = 2mc. The invariant mass m of j is m = rric iib — c and 
zero otherwise. In our approximation, X has zero invariant mass, M|- — {P — pY — 
The CM energy square of the viN collision is = {k + P)^. The four-momentum of 
W* is q = k — k'. As usual, we define = —q^ > 0, y = q ■ P/k ■ P, and the inelas- 
ticity variable z = pn ■ P/ q ■ P- In the rest frame of A^, y and z measure the relative 
lepton energy loss and the fraction of the W* energy transferred to H, respectively. The 
W*N CM energy square is — {q + Py — yS — Q^. The system X' consisting of 
j and X has invariant mass square M|^, = {q + P — pnY — ""^^ + ^ ~ z)S. 

As usual, we define the partonic Mandelstam variables as s = (g + p)^ = xyS — Q^, 
i = {q — phY = rn^ — xy{l — z)S, and it = {p — phY = M"^ — xyzS. By four-momentum 
conservation, we have s + i+u — M'^ + m'^ — Q'^. In the W*N CM frame, H has transverse 
momentum and rapidity 



Is (iu - MV) + g2 (tM2 + um? - 2M^m? 
s + Q' 

y*^-ln—. ^ ^ h-ln , (11) 

2 s (m^ - t) + Q\M^ - m^) 2s' 



respectively. Here and in the following, we denote the quantities referring to the W*N 
CM frame by an asterisk. The second term on the right-hand side of Eq. (11) originates 
from the Lorentz boost from the W*a CM frame to the W*N one. Here, y^ is taken to 
be positive in the direction of the W* three-momentum. 

We evaluate the cross sections of the partonic subprocesses (6)-(9) applying the 
covariant- projector method of Ref. [29]. In the case of process (6), our results can be 



5 



written in the form 



dy dQ"^ di 



-.{vi + a ^ r + cc[n] + 6) = 



y 



y 



+ (2 - y)AM 



(12) 



where FaWi ^aWi -^aW) and A^fn] are functions of s, £, u, and Q^, which are hsted 
in the Appendix. We combine the results proportional to the CO MEs iO'^ 



3p{8) 



and ( 



(8) 



exploiting the multiplicity relations of Eq. (5). By charge-conjugation 
invariance, the results for process (7) coincide with Eq. (12). The results for processes (8) 
and (9) may be obtained from Eq. (12) by scaling Fa[n] with the factor 1/2 to account for 
the spin average and by flipping the sign of as required by crossing symmetry. 

Invoking the factorization theorems of the QCD parton model and of NRQCD, the 
differential cross section of process (1) can be written as 



dydQ^ 



I- + H + j + X) = J^^ 
d^a 



dx 



dt 



a.b.n ^ ' dydQ^dt 



[Q2+(M+m)2]/(2/5) 

1^1 + r + cc[n\ + b) , 



(13) 



where fa/N{x, f^f) is the parton density function (PDF) of a in N, Hf is the factorization 
scale, (d^'a/dy dQ"^ di) (i/; + a — > Z~ + cc[n] + h) is given by Eq. (12), and 



t± 



1 

2i 



m\s - Q') -{s + Q')[s-M'^ JA(s, M^, m 



(14) 



with A(a;, y, z) — x"^ -\- y"^ + — 2{xy + yz + zx) being Kallen's function. The cross section 
formulas for processes (2)-(4) have the same form as Eq. (13). 

The kinematically allowed ranges of S, y, and are S > (M + m)^, (M + mY / S < 
y < I, and < < yS — (M + m)^, respectively. We then evaluate the cross section of 
prompt J/ip production as 

da{ui + N + ( J/V')prompt + j + X)^ da{ui + N ^l- + J/iP + j + X) 

2 

+ J2 B{XcJ J/^ + l)da{vi + N +XcJ + j + X) 
j=o 

+ B{ilj' ^ J/i/j + X)da{pi + N ^l- +1/;' +J+X). (15) 

The kinematic relations for uN and ep DIS given in Eqs. (15)-(24) of Ref. [18] and in 
Eqs. (17)-(24) of Ref. [16], respectively, carry over to the present analysis. In particular, 
they are not affected by a finite value of m. 
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3 Numerical results 



We are now in a position to present our numerical results for prompt J/il) production 
in CC DIS. To be specific, we consider i/Pb scattering at CHORUS and e~p scattering 
at HERA and THERA. The CHORUS experiment uses a wide-band neutrino beam that 
mainly consists of z/^ neutrinos, with an average energy of 27 GeV, and contains V^, v^, and 
Ve admixtures of about 6%, 1%, and below 1%, respectively. We evaluate the spectrum- 
averaged cross section per nucleon bound in lead using accurate parameterizations for 
the energy spectra of the individual neutrino species [22] as described in Eq. (38) of 
Ref. [18]. Furthermore, we approximate the effective nucleon PDFs of lead using Eq. (36) 
of Ref. [18], with A = 207.2 and Z = 82, and include a nuclear correction factor for 
incoherent scattering of 1.31 ^'20 [18]. At HERA, the lepton and proton energies in 
the laboratory frame are 27.5 GeV and 920 GeV, respectively. The respective values 
for THERA are planned to be 250 GeV and 1 TeV [19]. We adopt the same input 
parameters and PDFs, choose the renormalization and factorization scales, and estimate 
the theoretical uncertainties as in Ref. [18]. For simplicity, we take the CKM matrix to 
be diagonal. 



Table 1: Integrated cross sections (in fb) of prompt J/ip production by CC DIS in uPh 
collisions at CHORUS and in e~p coUisions at HERA and THERA. In each case, the 
NRQCD results for charmless and charmed X' and the CSM result (for charmed X') are 
compared. For comparison, also the respective LO results for fully inclusive CC DIS are 
presented. 





charmless 


charmed 


CSM 


inclusive 


CHORUS 
HERA 
THERA 


4.9+23 X 10"^ 
1 5+0.5 
■•-•^-0.4 
45+1.5 

^■'-> _12 


1.8 X 10-5 

0.50+S:- 

9 c+O.S 


1.7+;-^ X 10-5 
Q 43+0.14 


235 ± 39 
70.2 X 103 

195+^ X 10^ 

—5 



As a first step, we consider the integrated cross sections for CHORUS, HERA, and 
THERA, distinguishing between charmless and charmed underlying events X' and, in the 
latter case, between NRQCD and the CSM. The results arc presented in Table 1 together 
with their theoretical uncertainties. For reference, we also quote the respective LO values 
for fully inclusive CC DIS, evaluated from Eq. (14.8) of Ref. [28], choosing yu./ = 
where ^/ is an unphysical scale parameter of order unity, and otherwise adopting the 
theoretical input described above. We estimate the theoretical uncertainties on these 
values similarly as in Ref. [18]. We observe that, in the CHORUS, HERA, and THERA 
experiments, 27%, 25%, and 36%, respectively, of the prompt J/ip mesons are produced 
in association with an additional charmed hadron. The bulk of such events, namely, 94%, 
86%, and 80%, respectively, can be accounted for by the CSM. On the other hand, the 
CSM cannot explain, at LO, J/ip mesons associated with charmless underlying events. 
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Therefore, the observation of such events with the predicted rates would provide another 
piece of evidence for the failure of the CSM. As for CHORUS, the fraction of CC DIS 
events that contain a prompt J/ip meson is predicted to be 2.9 x 10^^ in NRQCD and 
7.2 X 10"^ in the CSM. Both values are compatible with the trimuon rates reported in 
Ref. [24]. As for HERA, the luminosity integrated over the entire operation period is 
expected to be of order 1 fb~^, so that there is not much phenomenological interest in a 
detailed study of the differential cross sections. 

In the following, we focus our attention on CHORUS and THERA. We choose the 
laboratory frame and take the J/ip rapidity yj/^ to be positive in the flight direction of 
the incoming lepton. The px and yj/^ distributions for CHORUS are displayed in Figs. 2 
and 3, and those for THERA in Figs. 4 and 5, respectively. In each figure, the NRQCD 
(solid line) and CSM (dashed line) predictions are compared. As explained above, the 
CSM contributions only arise from Feynman diagrams of the type shown in Fig. 1 (b). 

As for CHORUS, the NRQCD and CSM distributions exhibit maxima close to 
Pt = 1 GeV and rapidly fall off with increasing value of px because of limited phase 
space. The fall-off is significantly stronger in the CSM because the final state then always 
contains an additional heavy quark. The NRQCD and CSM yj/^ distributions have their 
maxima close to 3 and 3.5, respectively. For yj/^^3.5 GeV, they are similar in shape, 
while the CSM result significantly falls short of the NRQCD one as yj/^ is decreased. 

In the case of THERA, the NRQCD and CSM pt distributions have their maxima 
close to Pt = 5 GeV and smoothly fall off with increasing value of pr- The NRQCD and 
CSM yj/^ distributions have their maxima close to and 1, respectively. In the forward 
direction, they are similar in shape, while the CSM result significantly falls short of the 
NRQCD one as yj/^p is taken out to negative values. 



Table 2: Fractions of prompt J/ip mesons in uPh CC DIS at CHORUS that are produced 
directly and through the feed-down from XcJ and ■0' mesons. In each case, the NRQCD 
results for charmless, charmed, and unobserved X' and the CSM result (for charmed X') 
are compared. 





charmless 


charmed 


unobserved 


CSM 


direct 


43% 


76% 


52% 


77% 


via XcJ 


34% 


3% 


26% 


2% 


via ip' 


23% 


21% 


22% 


21% 



It is interesting to study which fractions of the prompt J/ip mesons are produced 
directly and by the various feed-down channels. This is done for CHORUS in Table 2. We 
observe that, within NRQCD, the values for the direct-production and XcJ-decay modes 
significantly differ depending on whether X' is charmless or charmed, while there is hardly 
any difference in the case of feed-down from ip'. As for charmed X', the NRQCD and 
CSM values are practically indistinguishable. However, if X' is unobserved, then NRQCD 
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and the CSM lead to significantly different results. 



4 Conclusions 

We studied prompt J/ip production in CC DIS at LO in the NRQCD factorization for- 
malism, providing the cross sections of all contributing partonic processes in analytic 
form. We presented theoretical predictions for vN fixed-target scattering at CHORUS 
and for e~p colhding-beam experiments at HERA and THERA, including conservative 
error estimates. 

There are two distinct classes of final states, one excluding and one including charmed 
hadrons in addition to the J /ip mesons. The former one results from pure CO processes, 
while the latter one may, to a large extent, already be explained within the CSM. Obser- 
vation of J/tl) mesons without any additional charmed hadrons would, therefore, provide 
evidence in favour of the CO mechanism. If the systems accompanying the J/'^ mesons 
are not observed, then there are two indicators of the CO mechanism: (i) the cross section 
is typically increased by factor of 4 compared to the CSM; and (ii) the fraction of J/ip 
mesons from XcJ decay is dramatically increased w.r.t. the CSM, while the one of direct 
J/ij) mesons is substantially reduced. 
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A Partonic cross sections 



In this Appendix, we present analytic expressions for the coefficients Fa[n], ^^[n], I/aM, 
and Aa[n\ appearing in Eq. (12). In order to compactify the expressions, it is useful to 
introduce the Lorentz invariants s = 2q-p, t = —2p-p', and u = —2q-p', which are related 



to the partonic Mandelstam variables by 



Q^t = t 



m 



and u = u + Q 



m 



respectively. Furthermore, g' = and mw are the gauge coupling and mass of 

the W boson, respectively, Gp is Fermi's constant, and Vu^ is the CKM matrix. 



vid{u) 
vid{u) 
vid{u) 



I CC 



(1)1 



CC 



CC 



3p(l)' 



6r>y 



3p(l)- 



u{d): 
u{d): 
u{d): 



F = T = L = A = 0. 



T 



(A.l) 
(A.2) 
(A.3) 
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vid{u) 
vid{u) 
vidiu) 



I cc 



I cc 



I cc 



3pf )] u{d): 
^Sf^\u{d): 
^Sf^]u{d): 



T 



T 



(A.4) 
(A.5) 



ud\ 



T = -2QH'^{2s + t) - 2Q^s[s^u - st(3t - 2u) - 2t^] - Q'^s\s^{t - 2u) + 2s{t^ 

-Atu - u'^) + t{2t^ - 2tu - u^)] - s^u{s'^ + 2st + 2t^ + 2tu + u^), 
L = 2QH{Q^-sf{s + t)\ 

A = ±{2Q^sH - 2Q^st{2s'^ + st - tu) + Q'^s^[s^{t - 2u) + 2s{f + u^) - tu{2t - u)] 
+s^u{s-u){s + 2t + u)}, (A.6) 



where the plus (minus) sign refers to the ui + d {ui + u) initial state. 



i'id{u) — > I cc 



3p(8) 



u{d): 



T 



(A.7) 



vis — > I cc 



c: 



277rMs4(g2 



m 



w) 



H2Q2 



u 



T = 32g^"(2s'' + UsH + 18sr - 39^) - 16g^[31s^ - s''(133t + 50m) + sH{328t + u) 

-3st2(33t - 32m) + 81t^{t + u)] + 2AQ^[A6s^ ~ s^{201t + 4m) + 28^(66^2 - 97tM - u^) 
-sH{10At^ + AOtu + 53^2) + Ust'^it + u){3t - 7u) - 15t^(t + uf] - 4Q^[244/ 
-4s^(173t - 57m) - 3/(61t2 + 168tM - 73m2) - s-'^(677t^ - 72t\ - 591tM^ - 154m3) 



+sH{t + M)(52t2 - 283tM + 259m=^) - 9st'^{t + ufiSt - 16m) + 3t%t 



+ 2g2s[152/ 



-52s^(5t - 2m) - 2/(131^2 - 16tM + 93m2) - s''^(469f'^ - 456^^^ - 531^^^ + 238m^) 
+4s2(t + u){25t^ + 7t\ + 56tu^ - 25u^) - st{t + uf{m'^ - 46iM+ 115m^) 
-6t^u{t + uf], 

L = -16g^°(4s3 - AlsH + 36si^ - m^) + 16g^[30s^ - 5^(73^ - 44m) + 5^^(55^ - 38m) 
-2,sf{m - 32m) + m^{t + m)] - 4g6[260s^ + s^(119i + 348m) + 28^(144^^ - 31iM 
-4m^) + ZsHillt^ + 166iM + 12m^) + 72sf{t + M)(3t - 7m) - 90t^(t + m)^] 
+4g^[220s^ + s^(391i + 150m) + 18s^(23i^ - IOm^) + s^(295i^ + I2>2t'^u + 18iM^ - 92m^) 
+sH{t + u){AQf - 2ntu + 175m2) - Qsf{t + uf{?>t - 16m) + 2>t^{t + m)^] - g2s[240s'' 
+4s^(197i - 40m) + 4s^(277i^ - 116iM - 183m^) + 3s^(293t^ - llOt^M - 143iM^ - 104m^) 
+2s2(i + M)(130i^ - 81i^M + 15iM^ + IOm^) - st{t + M)^(32i^ - 104iM + 269m^) 
-I2eu{t + m)^] - s^[lQ>s^ - As^{2m - 18m) - As^ilbt^ - 26tu - llu^) - s^(215t^ 
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-258t'^u - 2mu^ + lOAu^) -s'^{t + u){A7f - 22At'^u - 67tu'^ + 132u^) 
-su{t + uf{22t^ - mu + 40w^) - ?>tu^{t + u)\ 
A = -32g^°s(15s^ - 43si + 33i^) + IGQ^sflOSs^ - s^{m - 63ii) + st{A^t - Z9>u) - 9f(22t + 9u)] 
-24gS[90s^ - 5s^{t - 6u) + s'^{5t + 3u){37t - 5u) - st{155t'^ - UGtu - 76u'^) 
+21t2(t + u){3t + u)] + 4gS[300s^ - 8s^(10i + 9u) + 3s^(389t^ - UGtu - Ulu'^) 
-35^(40^^ - 215t\ + UOtu^ + 53u^) + st{t + u){253t'^ - 295tu - ITOm^) - 3t'^{t + uf{8t 
+5u)] - 2g^s^[120s^ - lOOs^t + 2s^(289t2 - 368tM - Qu"^) + 2,s^{77t^ + m'^u - 9btu^ 
+76ii^) - 2s{t + u){2t^ + 175^^-^ - 16^1*^ - Q2,v?) - t{t + ufi'it'^ - 28tu - 23u% (A.8) 



3p(l) 



c: 

12 



F 



^9'X\Vcs\ 



T = 256Q^6t(94s2 + 3Qst - 39t^) + 128Qi^[237s^ + s^{193t + 294m) + sH{58t + 637m) 
-3st^{39t + 8m) + 18t^(8t - 3m)] - 192Q^^[330s^ + /(668t - 57m) + s^(883t2 
+27UM - 426m2) + s2t(437t^ + 1039tM - 220m^) + 3st^{137t^ - 25tM + 30m2) 
+3t^(14t^ - llM^)] + 32Q^°[729s^ + s^(1805t - 3492m) + 6/(1234t^ - 1673tM 
-450m2) - 6s^(108t^ + 958t^M + 1907tM^ - 132m-^) + sh{63Ut^ - 5349t^M - 7635tM^ 
-1885m3) - 9st^(10U^ + 146t\ + 54tM2 + 30m^) - I2t\t + M)(38t=^ + 70tM - 13m2)] 
+ 16Q^[384s^ + s^(5636t + 813m) - s^(5617t^ - 22685tM + 6822m2) 
+6s^(978t^ + lieOt^M - 1498tM2 - 1763m3) - 2s3(410t^ - 13320t^M + 1257t2M2 
+12083tM^ + 1866m^) + sh{1568t^ + 4703f'^?7, + 1866t^M^ - 5U6fM^ - 4690m^) 
+3st^{t + u){339t^ + 86thi - 1166tu^ - 40m^) - 9t%t + uf{59f + 50tM 
+9m^)] + 24g^[840/ - 2s^(1954t - 1545m) - 5^(1024^^ + 3229tM - 3243m2) 
-s^(2423t^ - 536t2M - 7177tM2 + 2274m^) - s^(2897t^ + 4662t3M - 7860t2M2 
-1816tM^ + 5793m^) - s3(715t^ - 2492tS - 9804t3M=^ - 3932t'^u^ + 4957tu^ 
+2526n^^) - sh{t + M)(1832t'^ + 2228t^M - 2727tV - 39 Um^ + 1247m^) 
+3st^{t + M)2(286t3 + 92t^u - 365tu^ - 36m^) - 6t^{t + uf{64:t'^ - 39tu + Hm^)] 
-4Q^[4320s'' - 4s^(664t - 1305m) - 245^(454^^ - 107tM + 441m2) + 3s^(1112t=^ 
+6299t\ + 5326tu^ - 5709m^) - s^(14699i^ + 249i^M - 17283i^M^ - 6325iM^ 
-2196m^) - 3s^(2312i5 + 1739i^M + 3900i=^M^ + 11648i2M^ + 2576iM^ - 4419m^) 
+3s^{t + M)(1785i^ + 6260i^M + 5462t^M^ - 4586i^M^ - 2543iM^ + 1806m^) 
-sH{t + M)^(4141i^ + 1480i^M + 1815i^M^ - 1229iM^ - 1088m^) - 9st'^{t + uf{192t^ 
-182t\ + 51tu^ - 46m^) + 3t^{t + u)\l6t - 5u)^] + 2Q^s[336s^ + 40s^(184i - 39m) 
-4s^(1646t^ - 2065iM + 801m^) - 6s^(1326i^ + 2625t^M + 832iM^ - 571m^) + s^(4835i^ 
+18963i^M - 6609i^M^ - 10807iM^ + 7482m^) - 2s^(998i^ - 7945i^M - 13506i^M^ 
-8Ut\^ + 3488iM^ - 153m^) - 6s^{t + M)(473i^ + 1925i^M + 1679iV - 962t\^ - 858tu^ 
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+637u^) + 2s'^{t + uf{12St^ + 983t\ - 4764iV - 334^^1^^ + 3056iii^ 
-690u^) - st(t + uf{2m^ - 5728t^u - m'^u^ + limu^ + ^bu^) 
-&t^u{t + u)\m - hu)\ 
L = -128Q^^(12s^ - b9sH + 725^2 - 78^^) + 64Q^^[226s^ - 2s^(17i - l2Q>u) + sH{Am - bbu) 
+6st^(39i + 8u) - 36i^(8t - Zu)] - 32g^^[1040s^ + 10,s-^(137t + h9u) + 12s^(89i^ 
+21tM - 69^2) + ?>sH{lQmt^ - 652tM + 245w^) - 185^^(137^^ - 2htu + 30ii^) 
-18^2(14^^ - llu^)\ + 16Q^°[1122s^ + 2s^(2615t - 1518w) - 6s^(317t^ - 22htu + 826ii^) 
+12s^(789t^ - l2>Slt\ - lltu^ - m^) - sh{75A7t^ - 75A5t\ - 2559tu^ + 419ii^) 
+ 18st2(101t3 + Uet'^u + hAtu^ + 30m''^) + 2At^{t + u){m'^ + 70tM - 13^2)] 
+8Q^[3176s^ - s^(10037t - 12174m) + 2s^(8093t2 - 2803tM + 2982^2) - 2s^(10451i^ 
-19740t^M + ?,mu^ + 3146?/) + 2s-'^(10685t^ - 5025t'^?7, - 1929t^M^ + 1145tM^ 
-1602?/) - s2t(10778/ + 4106f'^?7, - mu'^u^ - AUWai^ - 2659m^) - Qsf{t + u){2,m^ 

+mt^u - iimtu^ - A<du^) + i8t^(t + uf{m'^ + sotw + 9u^)] 

-4Q^[8944s® - 8s^(1612t - 1673?x) + 3s^(9355t2 - 6847t?i - 3274^2) _ 2s^(12211t^ 
-23634t2?i + 3837t?x2 + 10282?i3) + 2s^(13886/ - 9175f\ - 92Aidt^u^ + 212tM^ 
-2017?i^) + 6s3(943t^ + 5888/m + 1572f'^?/ - 4856t\,^ - 1069tM^ + 378m^) 
-?,sH{t + ?i)(3608/ - 268f\ - 5828t^?/ + 1462t?r'' + 1785m^) 
+ 18st^(t + ?i)^(286t^ + 92t^u - ?,Qbtu^ - ?,&u^) - m\t + ufi&At^ - mu + IIm^)] 
+2g^[6912s^ - 8s*^(1121t - 408?x) + 6s^(2257t2 - 4629tM - 3094?/) - 3s^(655t^ 
-12272/m + 779^^2 + 3226?/) + 2s^(361l/ - 1259U^m - 4287t^M^ + 12773tM^ 
+8490m^) + 6s^(5158t^ + 5319/m - 7411tV - 7604tV + 2031t?i^ + 2315'u^) 
-6s^(t + ?i)(2059t^ + 892/?x + Sm^u^ - im'^u^ - 850tM^ - 358m^) 
+sH{t + m)2(6110/ + 7808t^M - 16386/m^ + 1526tM^ + 4139m^) - ISst^it + uf{im^ 
-l%2fu + hltv^ - A&u^) + Qt^{t + 'u)^(16t - buf] - Q'^s[992s^ - 32s^(55i + 39m) 
-45^(629^2 + 22bbtu + 762^2) + 4s^(1480t3 + lAlt^u - lUmu^ + 1873w^) + s^(173t^ 
-9t^u + 1887*^1*^ + 6293iii^ + 10272ii^) + 4s^(2459i^ - ^7m% - 5061i^ii^ + 2339*^-^^ 
-Amu^ - 1266w^) + 2s^{t + u){A9>22t^ - 92d>t^u - d>2m^u^ + l2Am'^u^ + imitu^ 
-hQMvP) - As^{t + uf{7m^ - 29ht% + 2937iV - ZQAQfu^ - ^mu^ 
+783w^) + st{t + 'u)^(3568i^ + m2t\ - 876fu^ + 13SAtu^ - 1391ii^) 
-12t'^u{t + u)\l6t - 5uf] - s'^[6As^ + Ss^{Alt + 38u) + 12s^(49i^ + 67tu + lOu^) 
+2s^(103i^ + 17Afu - 201tu'^ - 2A8u^) - s^(169i^ + 283t^u + 1995iV + lAAltu^ 
-152u'^) - 3s^(137i^ - 240^^-^ + 494iV + 896^^1*=^ - SBltu"^ - SOAu^) - 2s^(t + u)(820t^ 
-922t% + A9t^u^ + 1420^2-^=^ - lAAltu"^ + 172u^) - 2s^(t + u)^(6A0t^ 
-1288t'^u - 239tV + 1396fu^ - IbAOtu^ + 616u^) + 3su{t + 1^)^(464^^ - 216t^u 
-180t\'^ + 377tu^ - 168ii^) - 3tu^(t + u)\l6t - 5u)% 
A = 256g^S(35s^ + 107si + 141^^) - 128Q^^s[115s^ - s'^{t + 260u) + st{862t - 839u) 
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+81t'^{5t - 2u)] + 192g^2s[143s^ - 115s^(3i - u) + 5^(192^^ - 901tu + 359u^) 
-st{17St'^ + 1529tu - 269u^) - 3t'^{m^ + lOltu + 61u^)] - 32g^°s[2141s^ + s^(329i 
+3516u) - 6s^{255t^ - 537tu - 205u^) + 2s^{S26t^ + 765t'^u + 6279tu^ - 437ii^) 
-st{27Wt^ - SSt'^u - 7305tu'^ - 2645w^) + 3t'^{iS5t^ + 384i^w + 25Stu^ + 332u^)] 
+16Q^s[3304s^ + s^(7388i - 515u) + s^(559t^ + UOOdtu - IOTSSm^) + 2s^(3551i^ 
+1506i^w - 2277tu'^ - 5281m^) + 5^(466^^ + 14932t^'u - 108^^1^^ - 19618tu^ 
-3299ii^) + st(1783t^ + 2947^^-?/ + 362Athi^ - 3722tu^ - h&lhu^) - 27f^{t + u){At^ 
+87t\ + 27tu^ + 7u^)] + 24QS[440s^ - 2s^(1866t - 2009m) - 5^(2304^^ + mtu 
-4715m2) - s^(3041t3 + 2m'^u - A^mu^ + 2728'u^) - 8^(3933^^ - 430t^M - mm'^u^ 
+ 1020t?r'^ + 6377m^) - s^(351t^ - 3QUt*u - 9248t^M^ - 4764t\^ + 3517t?/ + 2512?/) 
-st{t + M)(2140t'^ + 1426t^M + 5UV - 989tM''^ + 1081?/) + 3t^{t + uf{At - ^u){2Qt^ 
+37tu - 10?/)] - 4Q/s[4336s^ - 4s^(332t - 1761m) - 24s^(488t^ - 207t?7, 
+383?/^) + s^(2936t^ + 14697/?/ + 7338t??2 - 19759??^) - s^(15136/ - 9795t^?/ 
-imAlt^u^ - 2225tu^ - 303??^) - 3s^(2535t^ + 2685/?/ + 4606t^??2 + 8496/??^ 
-57tu^ - 4637?/^) + s^{t + M)(2488t^ + 21065/?/ + 4390/?/^ - 13118/?/^ - 679Qtu^ 
+5749?/^) - st{t + ?/)^(4952/ - 1612t''^?/ - 933t2??^ - 901tu^ + 29u^) 
+3// + ?/)^(352/ - 528/?/ + 372tu^ - 71?/)] + 2Qls^[368s^ + 8/(956t - 179?/) 
-4s^(1510/ - 2033tu + 921?/^) - 2s^(3286/ + 7503/?/, + 3348t?/^ - 1085?/^) 
+5^(3713/ + 19387t''^// - 2925t^?/2 - 9199t?/ + 7240?/^^) - 4s^(1274/ - 1579/?/ 
-4845/?/^ - UbOt^u^ + 737tu^ - 267?/^) - 2s^{t + ?/)(932t^ + 2224/?/ + 7253/?/^ 
-1354/?/^ - 2087t?/^ + 1826?/^) + 2s(t + ?/)2(800f'^ + 3504/?/ - 3020t\^ - 322t^u^ 
+2448t?/^ - 767?/^) - t(t + ?/)^(784/ + 544/?/ - 1992^2^2 + i840t?/^ - 377?/^)]. (A.9) 



3d(1) 



UlS ^ I cc •^Pi 



c: 



8^?'^«2|v;,| 



277rM3s4(g2 + ^2^)2(2g2 + 5 _ 2t + M)4(2g2 _ 2s + t + ?/)5 ' 
T = 256g^^(2s=* + 18sH - 57s/ + 78/) - 128Q^^[4s^ + 5^(352^ + 13?/) - 3sH{53t - 60u) 
+3s/(108t - 19u) - 9/(16t + 29m)] - 192g^2[28s^ - s^(451i - 53m) + s'^(67/ + 152tM 
+97m^) - 3s^/26l/ + 104iM - 35m^) + 6s/(17/ + 19iM - 40m^) - 3/(28/ + 55tM 
+33m2)] + 32gi°[466s^ - 3s^(467i - 247m) + 3s^(209/ + 1163iM + 46m2) - s^(2735/ 
-4422/m - 3189iM^ + 488m^) - 27s^i(10/ - 153/m - IOO^m^ - 35m^) - 3s/(176/ 
+81/m - 381iM^ - 412m^) + 3/(i + m)(106/ + 161iM + 46m^)] - 16g^[929s^ 
+s^(1181i + 539m) + 3s^(735/ + 1277iM - 680m2) + s^(2573/ + 6228/m - 4521iM2 
-2440m=^) + 5=^(3920/ + 11969/m - 6441/m^ - 879Um^ - 628m^) - 3s^/148/ - 498/?^ 
+201/m^ + 1405iM^ + 900m^) + 3s/(i + m)(225/ + 55/m - 454iM^ - 293m^) 
-9/(i + m)^(22/ + 2Um + 8u^)] + 24g6[242s^ + 2s^(485i - 109m) + s^(1373/ - 467iM 
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-772u^) + s^{m9t^ - 1833t\ - 2817tu^ + 629u^) + s^(2965i^ + lU7t\ - 527UV 
+mtu^ + 1654ii^) + s^{107t^ - 1925t\ - 5218iV - 77St\^ + 2653iw^ + 713u^) 
+3sH{t + u){303t'^ + ISU^u - 250t\'^ + llStu^ + 262u^) - 3st^{t + uf{22f - 22tu 
-AW) + ^t^{t + uf{22t^ - 20tu + 9u^)] - 4Q^[76s'^ + 12s^(113t - A9u) + 3s^(879i^ 
-1226tM - M^) + s^(2963i^ - 8481^^-^ - lUStu'^ + 1745^3) + 3s^(695t^ - mit^u 
-UUfv!^ + 2915tu^ - 3A6u^) + 3s^(2320t^ + 493t^'u - 250^3-^^ + A396t\^ + UdOtu^ 
-1309u^) + s^{t + u){8m^ - 5620tSi - 566thj^ + 5962t\^ - 664t?/ - 1805?/) 
+9sH{t + ufiim"^ - 88t\ - 39t^u^ + 29tM^ - 85u^) + 3st'^{t + uf{4:f - 138t\ 
+llltu'^ - 62u^) - 33t^{t + u)^{2t - uf] - 2Q^s[88s^ + 8s^{17t + 23u) + 6s^{7t'^ 
+136tu - 7u^) - 4s^(265t^ - 588t^M + 237tM^ - 14m^) - s^(499t'^ - 2141f'^?i + 5307i2'u2 
+ 1393tM-'^ - 434?/) - 6s^(180t^ - 937/m + 94t^?/ + 768/m^ - 454tM^ + 61u^) 
-2s^{t + m)(1411/ - 169/m + 769tV - 296tV - 1819t?/ + 556m^) 
+2s'^{t + uf{A52t^ + 889t\ - 658t\^ - 178t\^ + 502tu^ - 253u^) - 3st{t + uf{32t'^ 
-120t\ + 150tV + 2tu^ + 15u^) - 66t^u{t + u)^{2t - u)% 
L = 128g^^(21s^ - 52sh + lUst^ - 156t^) - 64g^^[196/ - 2s3(142t - 31m) + 5s^t(179t 

-7m) - 6st2(108t - 19m) + 18t^(16t + 29m)] + 32Q^^[740s^ - 2s^(314t - 75m) + s-'^(2217/ 
-956tM - 401m^) + 6sh{76t^ + 68tM + 53m^) + 36st^(17t^ + 19tM - 40m^) - 18t^(28t^ 
+55tM + 33m^)] - 16g^°[1428s^ - 4s^(259t + 123m) + 6s^(58t^ - 689tn - 293m2) 
+6s^(586t^ - 609t'^u + 213tM2 + 42m''^) + s2t(317t^ - 2109t2M - 4809tM^ + 281m^) 
-6s/(176f^ + 81/m - 381tM2 - 412m''^) + 6t\t + m)(106/ + 16Um + 46m2)] 
+8g^[1429s^ - 2s^(526t + 827m) - 2s^(3841t2 + 2623tM + 855m^) - 2/(151t^ 
+6021t\ - 3543tM^ - 1217m-^) - s^(1718t^ + 7180t^M + 1764:t^u^ - 3364tM2 - 791m^) 
+2sH{im'^ + 107t^M + 1245t2M^ + 1166tM''^ - 925m^) + 6st^{t + M)(225t^ + 55t\ - AMtu^ 
-293u^) - I8t^{t + uf{22t^ + 21tM + 8u^)\ - 4g6[584s^ + 2s^(847t - 489m) 
-s^(13833i^ - 2389tM - 1102m^) - 8s^(1412t^ + 1143t^M - I827tu^ - 430m^) 
-2s^(3692i^ - 615t^M - 9bb2t^u'^ - 2788tu^ + 117m^) - 2s^(158U^ + 429U^m - 157U^m^ 
-3201i2M^ + 1801iM^ + 505m^) + 3sH{t + M)(10i^ - llGOi^M - 694^^-1^2 + QQ2tu^ + 789m^) 
-18st'^{t + M)^(22t^ - 22tu - 4lM^) + 18f{t + uf{22t'^ - 20tu + 9m^)] 
-2g^[72s'^ - 44s^(127i + 3m) + 3s^(2743i^ - 3472tM - 247m^) + 12s^(1199i^ + 26U^m 
-471iM^ + 139m3) + 2s5(573U^ - 9729i^M - 9186t\'^ + 546Um^ + 2391m^) + 6s^(1729i^ 
-92t^M - 4442t^M^ - 530t^M^ + 2035tM^ + 448m^) + 3s^{t + u){978t^ + 1200i^M - 1464i^M^ 
+60At\^ + 581iM^ + 37m^) - 2sH{t + M)^(205i^ + 352fu + 1428i^M^ - 548iM^ 
-695m^) - 6st'^(t + u)^(At^ - 138t\ + llltu^ - 62u^) + 66t^(t + u)'^(2t - u)^] + Qh[32s^ 
-56s^(83i - 5m) - 4s^(214t^ + 1993iM - 333m^) + 8s^(785t^ + 363i^M + 195iM^ + 209m^) 
+s^(8713i^ - 12949i^M + 593U^m^ + 816Um3 - 1288m^) + 4s^(2149i^ - 3863i^M - 4827i^M^ 
+3439t^M3 + 364iM^ - 1026m^) + 2s^(t + M)(1810i^ - 2224i^M - 1073iV + 657At'^u^ 
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-mtu'^ - imu^) + As'^it + uf{224:t^ - 207t\ - 539iV + 7m\^ - 2S2tu'^ - 158u^) 
+st(t + u)^(3AAt^ - 592t^u - 60t\'^ + 58Atu^ - A21u^) - 132t\(t + u)'^{2t - u)^] 
+3^1^^ + Ss^{83t + 9u) + 2s^(375t^ + 436tii - Uu^) - 16s^{7t^ + 69t\ + dltu^ + u^) 
-s^(117U^ - 129t\ + 407U2ii^ + A69tu^ - 78u^) - s^(1545t^ - 282At% - 93At\^ 
+Um'^u^ - 2179t«^ + 248'u^) - 2s^{t + u){A00t^ - 1450t^M + 79t\^ + 6A6t'^u^ - 1336tu^ 
+325m^) - 2s2(t + uf{9At^ - 518fSi + A89t^u'^ + 574:t\^ - Qlltu'^ + 240m^) 
+su{t + uf{80t^ - 184:thi + m\,^ + 185tu^ - I18u'^) - 33tu^{t + u)\2t - u)% 
A = -768g^S(s^ + lOst + 14^2) + 384g^^s[lls^ + s^(54i + u) - st{20t + 107u) + 3i^(60t + 7u)] 
-576g^^s[17s^ + 2s^(44t + 13m) - 2s^{Alt'^ + 36tu - Uu'^) + 2st{65t^ - 73tu + Ihu^) 
-t\76t^ + 202fAi + 69u^)] + 96Q^°s[112s^ + /(943t + 237m) - s^(187t2 - 102tM + 39^^) 
+s^(1213t^ + 552t\ + 957tu^ - 182u^) - st{782t^ + 71U\ - 1560tu^ - 481m^) + t^(236t^ 
+603fu + 7h9U? + 338?/)] - 48g^s[83s^ + s^(1527t - 47m) + 5^(1195^2 - 3Uu - 852m^) 
+s2(2387t^ + h7f^u - 1557tM2 - 835m^) + 5^(951^^ + 805t\ + 153t2M^ - 2199tM''^ - 167m^) 
+st(664t^ + 1522f^n + AQbfv^ - 1865tM^ - 1034m^) - 3t^{t + u){72t^ + ll8t^u + 83tM2 
+64m^)] - 72gS[26s^ - 2s^(127t - 103m) - s^(919t^ - 345tM - 268m2) - s^{<6AU^ - 943t^M 
-741tM2 + 271m3) - 2s^(540t^ - 440t^M - A9bt^u^ + 205tM^ + 299m^) - 5^(36^5 _ 279t^M 
-1294t^M2 - 482t^M''^ + 735tM^ + 239m^) - st{t + M)(62t^ - 12t^M - im'^u^ + 68tu^ + 237m^) 
-t'^{t + M)^(46t^ + 26t^M + 37tu^ + 3m''^)] + 12g^s[140s^ + 4s^(83t + 93m) 
-s^{U87t^ - 1290tM + 255m2) - s^(1337t^ - 1623t\, + 303tM2 + 959m^) - s^(1124t^ 
-7059t^M - 903t^M^ + 2507tM^ - 273m^) - s-^(282U^ + 663t^M + 900t^M^ + 38Ut^u^ 
+483tM^ - 1383m^) + s^{t + M)(224t^ + 2029t^M - 154t^M2 - 1702t2M^ + 247tM^ + 638m^) 
-st{t + M)^(484t^ - 320t''^M - 117t^M2 + 79tM''^ - 122m^) + t\t + M)'^(124t''^ - 84:t\ 
+72tu^ - 17m^)] - 6Qh^[A0s^ + 8/(39t + m) + 2/(151t^ + 208tM - IIIm^) - 45^(205^^ 
+12t\ + 81tu^ + 34m^) - s^(663i^ - 643i^M + ISSU^m^ - 129tM^ - 220m^) - 4s^(29i^ 
-757i^M - 57tV + 3m\^ - 298tM^ + 6m^) - 2s'^{t + M)(376t^ - n8t\ + 379tV 
+76iV - 448iM^ + 187m^) + 4s{t + uf{60t^ + 122t% - llSi^M^ + 21t\^ + 7ltu^ - 46m^) 
-t{t + ufiAQt"^ + 112i^M - 9m'^u^ + lOOiM^ - 29m^)]. (A.10) 



ViS I cc 



c: 



~ 1357rM3s4(g2 + rnlrf{2Q'^ + s - 2t + M)4(2g2 - 2s + t + m)^ ' 
T = 256g^^(30s^ + 2sH + 117si^ - 384i^) - 128g^^[306s^ - s^(608i + 87m) + s^i(2167i - 62m) 
-3sf{33&t - 337m) - 9t^{l88t - 57m)] + 192g^2[432s^ - s^(859i - 39m) + 5=^(421*^ 
-1574iM - 99m^) + sH{37hht^ + 1186iM - 223m^) - 12si^(151i^ - 200iM + 60m^) + 3t^{272t^ 
+729iM + 217m^)] - 32g^°[2784s^ - s^(1573i - 999m) - 3s^(6355i^ + 210Um - 102m^) 
+3s^(7219i^ - 8714i2M + 875iM^ + 292m^) - sH{<o87At^ + 6951i^M + 26382iM^ - 911m^) 
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+9st'^{AAAt^ + AlU^u - 2imtu^ - 460w^) + ?>t^{t + u){UU'^ - 1507iw - 1754^2)] 
+16Q^[2673s^ + s^{mm - 1239ii) - s^(44033i^ - I'^Ubtu + 41041^^) - 3s^(3193i^ 
+10588^2-^ - 12485iw^ + 22%u^) - s=^(27296i^ + 39099^=^1* - 21243iV - 26065^1*^ - Ih^u^) 
-s^i(3176i^ + 11324t=^u + 15915^^^^ - 37489iw^ - 4216w^) + 2,sf{t + u){lbhlt^ - 222ht^u 
+Umu^ + 4615ii^) - m^{t + uf{b2t^ + 2mu - b22u^)] + 24g^[126s^ - 2s^(2477i 
-1167w) + s^(12919i^ - 6887tM + 39421^^) + s^(9191t^ - 9%h7t^u - Zmitu^ + 2889ii^) 
+s^(16733t^ - 11163t^M - 34293^^1*2 + 9293^1*^ + 1848ii^) + s^(10549i^ - h2lt% 
-15282iV - llOQOt^M^ + 12535tM^ + 693m^) - sH{t + 'u)(1033i^ - 2^m\ + 7380^^1*^ 
-3908tM3 - 3512m^) + ?>st^{t + uf{<62t^ + 460^2^ - lOSltw^ + MISm^) 
+3t^(t + uf{im'^ - 312tM + 187^2)] - 4Q^[2868s^ - 4s^(349U - 1809m) + 3s^(3113t2 
-5398tM - 87?^^) + 3s^(3359t^ - SGGSt^w + ^Q^v"^ - 5429?/) + s^(22193t^ - 83613t^M 
-Ql'dUhP + 50309tM^ - 19170?/) + 3s^(19832t^ - llOll/w - 8238/^^ + 17540t2?7;'^ 
+6170tM^ - 3165m^) + 3s^(t + M)(1493t^ - 2368/m - 4062t^M2 + 14602t\^ - 3930tM^ 
-miu^) - sH{t + m)2(5225/ - 2788/m + 10263tV - 10177tM^ + 8483m^) 
+9s/(t + uf{im^ - 662t\ + 909tM2 - 486^3) - 3/(t + m)^(88/ - 136tM + 55u^)] 
+2Q^s[1656s^ - 8s^(767t - 381m) - 2s^(2915/ + 5108tM + 297^^) + 12s6(305t^ 
+873t\ - 409tM2 - 105m^) + s^(5623/ + 2463/m + 43455t2?/ - 13547tM-^ + 6006m^) 
+2/(9500/ - 26299/m + 6042/^^ + U996t\^ - 13858tM'^ + 4347m^) + 6s%t + M)(1585t^ 
-3807/m - 2069/^2 + 3540t\-'^ - 3635t?/ + 666?/) - 2s^{t + m)2(3368/ 
-6583/m + 1098/^2 - 154/m^ + 818^??'^ - 291m^) + st{t + m)''^(1016/ - 4048/m 
+719 - 5674tM^ + 2345m^) + 6t\{t + u)\88t^ - 136tu + 55u^)], 
L = -128g^^(9/ - 286s2t + 234s/ - 768/) + 64Q^^[302s^ - 2/(892t - 201?/ - sH{329t 
+653m) - 6s/(336t - 337m) - 18/(188t - 57m)] - 32g^2[2164/ - 26s^(61t - 104m) 
-3s=^(1565/ - 456iM - SIm^) + 65^^(304/ + 917tM + 140m2) - 72s/(15l/ - 200iM 
+60m2) + 18t3(272/ + 729tM + 217m2)] + 16g^°[6378s^ + 4s^(3266t + 1281m) 
-6s^(208/ - 2373iM + 671m2) + 6s^(3044/ + 965/m + 877tu^ - A38u^) - 5^^(4111/ 
-34593/m + 18147iM2 - 4367m^) + 18s/(444/ + 475/m - 2109tu^ - AGOu^) 
+6/(t + m)(346/ - 1507tM - 1754m2)] - 8g^[7745s^ + 4s^(10507i - 726m) + 2s^(16919/ 
+7391iM - 10689m2) + 2s^(27235/ - 7197/m - 8715tu'^ - 4997m^) + s^(35626/ 
+25776/M - 17340/m2 + 40iM^ + 627m^) + 2s2i(5453/ - 7960/m - 31227/m2 
+14432iM^ - 1234m^) + 6s/(i + m)(1551/ - 2225/m + 3494iM2 + 4615m^) 
-18/(i + m)2(52/ + 277iM - 522m2)] + 4g^[88s^ + 2s^(23351i - 10201m) + 3s^(25513/ 
-6487iM - 9136m2) + 45^(23401/ - 21672/m - 1643Um2 + 1690m^) + 2s^(4768l/ 
-16807/m - 15504/m2 - 6196tM^ + 9008m^) + 6s^(8257/ - 3337/m - 4869/m2 + 5617/m^ 
-3265iM^ + 723m^) + 3sH{t + m)(6034/ - 10736/m + 10784/m2 - 5524iM^ - 3909m^) 
-18s/(i + m)2(62/ + 460/m - lOSliM^ + 1413m^) - 18/(i + m)^(134/ - 312tM + 187m2)] 
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+2g^[7656s^ - 4s^(274U - 5403ii) - 3s^(21617i^ - Umtu - IQlu^) - 6s^(15638i^ 
-19355t\ - UOStu^ + 6329w^) - 2s^(55853i^ - 72369^3^ - 33240i^w^ + 19817iw^ 
+133591^^) - 6s^(14998t^ - 13515t^ii - 6616tV + 6070^^-^^ + 4212iii^ - 1211^^) 
-3s^{t + 'u)(10778t^ - 24948i^w + 5428iV - 12752i^w^ + 8563iw^ - 945ii^) 
+2sH{t + ufimm'^ - 856t\ - 9375t\'^ + 6659tu^ - 7720m^) + 18st^{t + uf{im^ 
-662t'^u + QOdtu^ - mu^) - 6t^{t + uf{88t'^ - 136tu + 55u^)] - Q^s[4:672s^ 
+8s^(1867t + 675^) - 4s^(2290t2 - 6581tu + 2793it^) - 8s^(6407t^ - 5820^^-^ + 2055iw^ 
+1442ii^) - s^(72785i^ - 122829t3'u + 7539^2^^ + A7513tu^ - 17400w^) - 4s^(16577t^ 
-32131t^M - 2523t^'u2 + 22739^2^3 + 4412tM^ - 6702m^) - 2s^{t + M)(14446i^ - AOdUt'^u 
+33235t^M2 + 98m\^ + 2753tM^ - 5AUu^) - As^{t + uf{76t^ - I283t^u 
-ISm^u'^ + 1924t2?/ - 1754t?/ - 264m^) + st{t + M)^(2504t^ - 9568t^u + U628t'^u^ 
~11872tu^ + 5843?/) + 12t^u{t + u)^{88f - 136tu + 55u^)] - s^[8s^ - 8s^(509t - 79m) 
-6/(971t2 + 944tM - 259^2) + 4/(1348t=^ + 21t\ + 7mu^ + ISSm^) + s^(14551t^ 
-2921^2^ + 22191t\2 + 973tM^ - 602m^) + 3s^(3143t^ - 9528t^M + Um^v^ + ?,8A8fu^ 
-3933tM^ + 560m^) + 2s^{t + u){m2t' - 12622t^M + 7303t\^ + 7Qt^u^ - Q7mu^ + 1855'u^) 
+2s2(t + uf{m^ - 7A8t\, + 3199t^?/ + IQm'^u^ - 2971tM^ + 1138m^) 
+2>su{t + uf{<6m^ - lOSOt^M + 828^2^^ - 287tM^ + 150m^) + 2,tu^{t + uf{88t^ 
-imu + hW)], 

-256gi^s(47s2 - imst + 168t2) + 128Qi'*s[481s-'^ - s2(394t - 223m) - st(608t + 329m) 
-297^2^] - 192gi2s[623/ + 2s^(255t + 92m) - 2s'^{lA7t^ - 379tu + 35u'^) - 2st{6Alt'^ 
+46tM - 280m2) - 3t2(922t2 + 512tM + 121m2)] + 32Q^°s[3206s^ + /(8171t - 1281m) 
+3s3(1611t2 + 678tM - 1441?/) + s2(1636U^ + 22752/^ + 2829tM2 - 1294m^) 
-st(30490t^ - 16863/m - 11742tu2 - 623m-'^) + 3t^(332/ + 6279/m + 5109tM^ + 1322m^)] 
-16g^s[1447s^ + lls^(144U - 547m) + 5^(5425*^ - 10841tM - 7104m=^) + s='(72799t=* 
-4863t^M - 25425tM^ + 721m3) - 5^(5911^^ - 62125t^M + 37611iV + 19939iM=^ - 1171m^) 
-si(10766i^ - 37108i^M + 19803tV + 29429iM^ + 7112m^) + 27t'^(t + u)(16At^ + 56t\ 
-mtu^ - 182u^)] - 24g''s[866s^ - 2s^(2379i - 1481m) + s^(1161i^ + 1727iM - 142m2) 
-s^(22889t^ - 28777t^M + 595iM^ + 5245m^) - 2s^(10917t^ - 3350t^M - 14283t^M^ + 5145iM^ 
+1840m^) + s^(3411i^ - 8839i^M + 33722iV - 906f^M^ - 8671iM^ - 673m^) - st{t + M)(2542i^ 
+2224i=^M - 5790i2M^ + 6874iM=^ + 3031m^) - 3t^{t + ufimt^ - 1050t\ + 399tu'^ 
+307m=^)] + 4g^s[3388s^ - 4s^(497t - 1419m) + 3s^(3527t^ + 54tM - 217m^) - s^(17365i=^ 
-38457i^M + 23991iM^ - 1835m^) - s^(58828i^ - 103125i^M + 42771i^M^ + 32773iM^ 
-14379m^) - 3s=^(8829i^ - 1689i^M - 27392iV + 30108i2M=^ - 417iM^ - 4199m^) + 
s'^{t + u){388t^ + 2699i^M + 10138iV - 30338i^M^ + 9539iM^ + 2992m^) - st{t + M)^(6512i^ 
-21184i=^M + 16905i^M^ - 3787iM^ - 4648m^) + 3t^{t + uf{532t^ - 600i^M - 96iM^ + 199m^)] 
-2Q^s^[106As^ + 8s^(419i + 7m) + 25^(4483^^ + 2524iM - 2607m2) - 45^(1025^=^ + 3783i^M 
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-3645^1^^ + 1589w^) - s^(18649i^ - U3S9t^u + TlQgii^-u^ - 36imu^ + 4220ii^) 
-4s^(2465i^ - 12601^^-^ - 551UV + 13253t'^u^ - m2tu' + 1290ii^) - 2s'^{t + 'u)(1988i^ 
-4358i^ii - 5875^^-^^ + 3644^2-^^ - mA2tu^ + 2063ii^) + 4s(i + 1^)2 (136^^ 
+942t^'u - 4783^^1^^ + 2QAlt^u^ + mu^ - 259ii^) + i(i + uf{MAt'^ - n92t\ + lM%tW 
+620tu^- 853ii^)]. (A.ll) 



vis ^ I cc ^Sq 



lo(8) 



c: 



967rMs4(Q2 + m2^)2(2g2 ^ g _ 2t + uy{2Q^ -2s + t + u)^' 
T = 32gi°(2s3 + lOs^t - Gst^ - 3t^) + 16g^[5s^ - s^{15t - 26m) - s2t(44t + 17m) 
-3st=^(3t + 8v) - 9t^{t + u)] - 24:Q^[Us^ + s^{m + 16m) + 2s^{2t^ + 27tM + m^) 
+sh{2t + 5uf + 12st'^u{t + u) + 3t%t + uf] + 4Q^[44s® + 8s^(17t - 6m) 
+3s^{A9f + 20tM - 49m2) + ^^(Sgt^ + 12t^u - 135tu^ - 82u^) + sH{t + u){AQt'^ + 35tu 
-59m2) + 3st^{t + uf{t - 8m) - 3t^{t + uf] + 2g=^s[8s^ - 4s^(27t - 20m) 
-2s^(43t2 - 44t?7, - 15m2) - s^(6U^ - m^i - 3U? + 94m^) - 4s^(t + M)(3t^ - 17t\ 
-16tu^ + 13?/) - st{t + uf{2t'^ - Utu + llM^) - 6t'^u{t + u)% 
L = 16g^°(4s'^ - sH + 12st2 + 6t^) - 16g®[10s^ - s%15t + 8m) + sh{m - 10m) 
-3st^(3t + 8m) - 9t^{t + m)] + 12g^[8s^ - s^(4U + 20m) - 10s'\(2t + 3m) 
+sH{f + 26tM + 36m2) + 24st^M(t + m) + 6t^{t + uf] + 4g^[8s6 + s^{85t + 24m) 
+6/(2t + M)(llt - 2m) + s^{5t^ - A2t\ - IStu^ + 8m^) - sh{t + u){2t^ - 41tM - 79m2) 
-3st^(t + uf{t - 8m) + 3t^{t + uf] - g2s[32s^ + 4s^(21t - 4m) + 4s^(6U^ 
-8tM - 15m2) + s^{2l5t^ + ISOt^M - 69tM^ - 88m-'^) + 2s^{t + uf{m'^ - Wtu - 38u^) 
+st{t + uf{AOt^ + 8tu - 77m^) - I2t\{t + uf] + 3s^[4s^t + 4/(3t^ - m^) 
+sH{21t'^ - 6tu - 23m^) + s'^{t + M)^(13t^ - 29tM + 12m2) - su{t + uf{6t^ + tu- 8m^) 
+tu\t + u)% 

A = 32gi°s(s2 + list - 9f) + 16g^s[13s^ - s^(18t - 31m) - st{25t + 22m) - 9t^{2t + u)] 

-24gS[22s^ + s^{19t + 26m) + s^{17t'^ + A2tu + v?) + hst{t'^ + %tu + 4m^) + - u^)] 

+4gS[76s5 + 4s^(4U + 3m) + 35^(57^^ - Uu - b3u^) + 8^(100*^ + m\ - 23Atu^ - 95u^) 
+st{t + u){17f - mu - IOm^) + 9t'^u{t + uf] - 2Q^s^[8s^ + 4s^(33i - 10m) + 25^(41^^ 
-76tM - 21u^) + s'^{m^ - 126t\ - Ultu^ + 68m^) + 2s{t + u){6t^ - llfu + ?,2tu^ + 31m^) 
+t{t + uf{At^-Atu-nu% (A.12) 



vis l-cc fef ^ 



c: 



4327rM3s4(g2 + rnlrY{2Q'^ + s - 2t + M)2(2g2 - 2s + i + uY{-AQ'^ + s + i + 4m) 
T = -256g^^(68s^ - lOOs^t + 81st^ - 129t^) + 128g^^[754s^ - s=^(262t - 337m) 
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+5^^(1303^ - 323u) - 3st'^{2m - 7iu) - 18t^{8t + 43ti)] - 192Q^2[1058s^ + s^(201i 
+8Uu) + 5^(1686^2 - 12mu + 134x^2) + 5^^(485^^ + 2350tu + mu^) - 9st^{2St^ 
+n9tu - llSii^) + 3t^{31t'^ + 76tu + 33u^)] + 32g^°[6323s^ + 26s^(121i + 309u) 
+3s^(2281t^ - 5643tu + 482^^) + s^(5849f^ + 12807i\ - UUdtu^ - 3397u^) 
-sH{15Mt^ + U97t\ - 21309^1*^ + lUlu^) + 18st'^{2t + u){18t'^ + mu - 17u'^) 
+3t^{t + u){77f + 736tu + SUu"^)] - 16g^[5698s^ + 2s^(2839t + 7232ii) 
+s^(2665i^ - 30689iw + 4:128u'^) + /(1819t^ - 4218i^M - 46863iM^ - 5648n^) 
-2s^(3065t^ + 7mt^u - 5454i^'u2 + 17027tM^ + 2044m^) - 45^^(758^^ + 2174^^^ 
+3540t\2 + 241 ItM^ + 800m^) + 3st'^{t + u){270t^ + 2297t^u - 1223tu'^ - 1783m^) 
+9t^{t + ufiAt'^ - 174tM + lOlM^)] + 24Q^[364s^ + s^(2043t + 4031?;,) + s^{126t'^ 
-eemu + 3871^2) - 4s^(1090t^ + 3820t2M + 6803tu'^ - 28u^) - s^(116U^ + 5383t^'u 
-4443t^?7,^ + 10389t?i^ - 2834m^) - s3(4647t^ + 15326t^M + 18464f\2 + 6726t2M2 
-2320tM'^ - 2197m^) + sh{t + M)(37t^ + 973t^M - 7686^^^^ - 5822tu^ + 505m^) 
+9st^{t + uf{4:lt^ - 128t^M + 129tu^ - 987?) - 3t^{t + uf{m'^ + lOtu + 162?/^)] 
+4Q^[1288s^ - 4s*^(827t + 309m) - 24s^(28t^ + 78tM + SSIm^) + 2s6(4735t^ + 16281t\ 
+20409tu^ - 9989^3) + s^(9605t^ + 51879f\ + 52413t^M2 + 48737tM3 - 8922m^) 
+6s^(58U^ - 513t^M - 6977tV - 11518t2M^ - 6696t?/ - 1421m^) + s^{t + u){7366t^ 
+27743t^M + 47203tV - 10811t^M^ - 4083Um^ - 8162m^) 

-sH{t + M)2(3773t^ + 1034f\ + 4305tV - 3214tM-'^ + 6752m^) + 9sf{t + uf{UAt^ 
+28t^u - 6tu'^ - 571m^) - 3t%t + u)'^{38t^ - 20tu + Uu'^)] - 2Q^s[512s^ - 4/(113t 
-539m) - 4s^(124t2 + 1255tu + 1047m^) + s^(1267t^ + 5985t^u + 2457tu^ - 15461m^) 
+s^(4181t^ + 21109t^M + 22725t^M2 + 7343tM^ - 12422?v;^) + 2s'^(1088t^ + 4721t^'u 
+5460t^'u2 + 2156^2^^ + 3338tM^ + 2577m^) - 2s^{t + M)(203t^ + 4561fV + 11810^^^ 
+20702t'^u^ + 3997tit^ - 8011m^) + s\t + uf{772t^ + imm^u + 10004*=^^^ 
-UAQfu^ - 9h79tu^ + 8075m^) -st{t + uf{28t'^ + 2Am^u + IhhAfu^ + 3Aiotu^ 
-4901m^) + Qt^u{t + uf{38t^ - 20tu + 41^^)], 
L = -128g^^(41s=^ - 735^^ - 162sf + 258t^) + 64gi^[321s^ + s^{53t + 311m) - sH{722t 
-1243m) + 6si^(243i - 74m) + 36t^{8t + 43m)] - 32Q^2[841s^ + 7s^(388i + 189m) 
-s^(846i^ - 8636iM + 991m^) + 3s^i(1331i^ + 49iM + 853m^) + 545^^(23^^ + 119iM - 115m^) 
-18t^{31t^ + 76iM + 33m^)] + 16g^°[245s^ + 14s5(628t - 129m) - 18s\26t'^ - 1213tu 
+531m^) + s^{5597t^ - 104:lfu + 5769tM^ - 5131m^) + sH{4:562t^ + 15765fu - 26229iM^ 
+5273m^) - 36st'^{2t + u){18f + 43iM - 17m^) - 6t^{t + u){77f + 736tu + 344m^)] 
+8g^[3348s^ - s^(9239i - 19819m) + 25^(1279^^ - 9989iM + 18843m2) - 6s^(144i=^ 
-3613i^M - 382tM^ - 3980m^) - s^(6371i^ + 11822i^M - 75102i2M^ + 4946tM^ - 2395m^) 
+sH{2060t^ + 10469i^M + 6942i2M^ + 5018tM^ + 3137m^) + 6st'^{t + u){270t^ 
+2297t^u - I223tu^ - 1783m=^) + I8t^{t + ufiAt^ - l7Atu + IOIm^)] - 4g^[8312s^ 
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+s'^{im + 45201ii) + 5^(5844^2 + Smtu + Ibbllu^) + 4s^(2702i^ + 14361^^1* 
+3267^1*2 + m%2u^) + s\Amit^ + 30858i=^ii + Ull2m'^u^ - bl22tu^ - imbu^) 
+s^(3960i^ + 33346^^ + 54082^^-^^ + 42090tV + 7969iw^ - 8561w^) + ZsH{t + 'u)(1646i^ 
+9812i^ii - mAltW - 5333tu^ + 4713w^) + bAsf{t + uf{Alt^ - I2%t^u 
+129tii^ - 98m^) - m^{t + uf{A?,t^ + mu + 162m2)] + 2g^[9520s'' + 2s^(2231t 
+24279m) + 185^(335^2 + ISOltw + 43821^^) + s^(20915t^ + 73623^^^ + U^mu^ 
+25049^^) + s^(14459t^ + 66834^^-^ + UMA2tW - 51646^1^^ - 46257ii^) + 9s^(424i^ 
+5105t^'u + 9228^^7/2 + AZmt^v? - 3902tu^ - AAlQu^) + s^{t + 'u)(12028i^ + AQZ%2t% 
+2989tV - 33221tV + 15791tM^ - 8303m^) + sH{t + ufibQm^ 
+ 1756f\ + U'im'^u^ - 5237tu^ + 3391m^) - 18st'^{t + M)^(144t^ + 28t\ - 6tu'^ 
-571m^) + 6t%t + u)\38t'^ - 20tu + 41m2)] - gls[5424s^ + 4s^(307t + 6715?i) 
+4s^(818t2 + 3179tM + 10041n2) + 3s^(5167f'^ + 15237t\, + 13fM^ - 2217?/) 
+s^(8507t^ + 55445t\ + 87069t'^u^ - 68617tM^ - 64252m^) + 2/(2522t^ + 16451t^M 
+24828^3^2 ^ i754t\3 _ 28708t?/ - 23019m^) + 6s^{t + M)(1591t^ + 5269^^^ - 1673^3^2 
-9487t\^ + 1194tM^ - 592m^) + s^{t + uf{3976t^ + 11984t^M + 7A92t\^ 
-27182t\,3 - 15277t?/ + 3613m^) - st{t + ufilOOOt^ - 1136t^M - 43Ut'^u^ - 3362tu^ 
+11911?/) - 12tS/t + u)\38t'^ - 20fM + 41?/)] + s^[5Us^ + 4s®(43t + 831m) + 4:s\3t^ 
-mu + llUu^) + s^(1751f^ + 4593/m - 4467tM2 - 5821m-'^) + s^(1490t^ + 10845/m 
+ 15927t2M2 - 13873t?r'^ - 17133m^) + s^(741t^ + 3982t^?i + 8470t^M2 + 4356/m^ - 7895tM^ 
-7886m^) + 2s^{t + m)(590/ + 1567/m - 4234/^^ - 8668t2M''^ + 3472tM^ + 4391m^) 
+s'^{t + tif{562t^ + 2464/m - 89At\^ - 13070/m^ + 589tM^ + 9567m^) + su{t + m)^(524/ 
+848t^M - 1230t'^u^ - 1630tM^ + 2543m^) + 3tu^{t + uY{38t^ - 20tu + 4:lu% 
A = 256Q^^s{72s^ + 59st + 2A9t^) - 128g^^s[678s^ - s'^{32t - 69u) + st(1333t - 983m) 
-9t'^{5t - 76m)] + 192gi2^[909s^ - s^(277i - 594m) + 52(1018*^ - 2285iM + 138m2) 
-st{6mt'^ - 346tM + 889m2) - 9t'^{5Qt'^ + Ultu + 13u^)] - 32gi°s[5967s^ + s^(347i 
+9459m) + 6s%884t^ - 2756tu + 927u'^) - 3s^{3A9t^ - 3067i2M + 1579tu^ - 933u^) 
-st{5674:f + 12369i2M - 2715tu'^ + 337u^) - 3t\61t^ + 1299i2M + 2334:tu^ + 988m^)] 
+16g^s[7302s^ + 2s^(3653i + 10740m) + 5^(7738^^ - 14657iM + 20043m2) + 6s^(1062i=^ 
+2918i2M - 2711iM2 + 890m^) - 25^(1498^^ - 6786i^M - 305m\'^ - 9941iM^ - 1884m^) 
+si(247i^ - 9788f=^M - 10128iV + IISOO^m^ + 9325m^) + 9t^{t + u){83t^ + i02t\ + 231tu'^ 
+263m^)] - 24gS[1476s^ + s^(4237t + 8073m) + 5^(2864^^ + 1595tM + 10167m2) 
+s^(4117i^ + 6737i2M - 7049iM2 + 3303m^) + 2s^(170i^ + SlOOi^M + 9969i2M2 - 3700iM^ 
-2043m^) + 5^(2531^^ + 8522i^M + 20322iV + 22976i2M^ + 3530iM^ - 3171m^) 
+st{t + M)(463i^ + 1795i^M - 2562fu^ + 928tu^ - 443m^) + 9t'^{t + uf{10t^ 
+101t\ + 277tu^ + 24m^)] + 4gS[792s^ + 4s^(2045i + 3123m) + 125^(699^^ 
+1642iM + 1953m2) + 6s^(425i^ + 289t\ - 3805tu'^ + 927m^) + s^(8161t^ + 36723i^M 
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+38U5t\^ - 37163tu^ - 18666w^) + 3s^{629t^ + 3358i^ii + 7056fu'^ - 736t\^ 
-13541^1^^ - 6318u^) + 3s^{t + u){2236t^ + S557t\ + 3573t^u^ - 17i5t\^ - 9575iw^ 
-2154ii^) - st{t + uYimSt"^ + 1228t=^ii - 6141^2-^^ + 12280iii^ + 5U9u^) 
+3t'^{t + uf{104:t^ + 186t\ - 426iii^ - 211w^)] + 2Q'^s'^[96s^ - 4s^(355i + 303ii) 
-4s^(934t2 + 27S5tu + 2367u'^) + 3s^{171t^ - 623t\ + 1265tu'^ - 581u^) - s\569t^ 
-5535t\ - Xmm^v?- - 40237tM^ - 19446w^) - 4s^(890i^ + 5930t^M + 9786^^^^ 
+2525t^u^ - 3%77tu' - 1008xi^) - ^s^{t + u){7U^ - 261t^it - 2086i\^ - lYIU^u^ 
+1437tM^ + 2975'u^) - s{t + uf{2hU^ + 3316t^'u + 2856tV - 19172^21^2 - IZmtu" 
+8853ii^) - t[t + uf[mt'^ - 2m^u - 9m^u^ + 1294^1^2 + 1081ii^)]. (A.13) 



VlS I cc 



F 
T 



c: 



' 247rM3s4(Q2 + m^)2(2Q2 + ^ _ 2t + uy{2Q^ -2s + t + uf' 

-- -768Q^^{8s^ + 2sh - 16st^ - 29t^) + 384gi^[97s^ + lls2(llt + 2u) + sh{382t + 177m) 
-st^{21t - 184?7,) - 6t%76t + 11m)] - 576g^2[114/ + s^(184t - 33m) + s^{U7t^ - 421tu 
-134^2) + sH{U05t^ + 587tu - lUu^) - st\425t^ - 449tn - lOu^) + t^{242t^ + 596tu 
+207^2)] + 96g^°[389s^ + 5s^(157t - 132m) - 2s^(1916t2 + 2421tM + 162^2) + 28^(3282^^ 
-4698t2M - 1407^2 + 304^3) - sH{6m^ + 4893t\ + 8511^^ + 637^^) + st\ll39t^ 
+966t\ - 503Atu^ - 177Au^) - At^{t + u){8t^ + 355hi + 311m=^)] + 48Q^[68s^ 
-s^(136t - 725m) + s%10631t^ + 2109tM - 1206^2) + 2s^(2248t2 + 5992t2M - 6168^^2 
-1325^3) + 2s2(4294t^ + 9814:t^u - Amfu^ - 7imtv^ - 434?/) + sH{h8At^ + 4035^2^ 
+3342t\2 - 11746tM2 - 3258m^) - st^{t + M)(471f^ - l&22fu + 3626^^^ + 3676m-'^) 
-3t^{t + uf{29t^ - 126tM + 367^2)] - 72Q6[96s^ - 10s^(2t - 3m) + s^{39mt^ - lA83tu 
+u^) + s^(3451t2 - 3472t\ - 3477tM2 + 1314m2) + s^(6017t^ - 1198t2M - 12008t2M2 
+2188tM2 + 2157m^) + s''^(2495t^ - 1676/m - 79m^u^ - 3A2Qfu^ + 5105tM^ + 910m^) 
+sH{t + M)(580t^ + 1156f\ - 2443tV + lOGltM''^ + 1443m^) - st^{t + uf{98t^ 
-276t\ + i53tu'^ - 1036^2) + 2t^{t + M)''^(88t2 - 137iM + 75m2)] + 12g^[208s^ 
-4s^(476t - 159m) + 45^(126^^ - 1584iM + 279m2) + s%38AAt^ - 11529t2M + 2658tu^ 
-553m2) + s^(7955i^ - 27135i2M - 55231^"^ + 16307iM2 - 5196m^) + /(18632i^ - 5775i^M 
-4692iV + 2A368t'^u^ + 7272tu^ - 6201m^) + s^{t + M)(965i^ - 7412i^M - A9m^u'^ 
+15238iV - 2215tM^ - 2246m^) + sH{t + uf{l3t'^ + 2h&t\ - 2313tV + 2299tu^ - 2516m^) 
+sf{t + uf{h92t^ - I782t^u + 2mitu^ - lU2u^) - t^{t + uy{232t'^ - 232tu + 67u'^)] 
-6Q^s[272s^ - 8s\28At - 91m) - 45^(148^^ + 933iM - GSm^) + 2s^(1202t=^ + 1255i2M 
+A8tu'^ - 349m2) + s^(1045t^ - 2987t2M + 13929t2M2 - 881tu^ + 2Uu^) + 2s^(2582t^ 
-8601i^M + 30t^u'^ + 4:778t'^u^ - 3600tu^ + 1071m^) + 2s^{t + u){2155f - 2009i^M - 307t^u'^ 
+18A2t\^ - 3922iM^ + 1027m^) - 2s'^(t + u)^(106At^ - 76U^m - 724iV 
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-mt\^ + 856tu'^ - 310u^) + st{t + uf{320t^ - 1776t^u + lS90t\^ - 1066tu^ 
+555w^) + 2t\{t + u)^{232t'^ - 232tu + 67u% 
L = -128Q^^(40s=^ + 67sH + 96st'^ + 17At^) + 64g^^[266s^ + 10s^(6U - 6u) + sH(2153t 
+347w) + 6st'^{21t - ISiu) + 36t^{76t + llu)] - 32g^2[384s^ + 2s^(129t - 653w) 
+20s^(396i^ - 137^1* - 621^^) + 3sH{235t'^ - UGStu + lOTw^) + 18st^(425t^ - U9tu 
-10«2) - 18^^(242^2 + 5mu + 207m2)] - 16Q^°[1022s^ + 2s^(4519i + 2706m) 
-6s^(149i2 - 28S5tu - iUu'^) + 4s^(497t^ + 4173t\ + 2mu^ - 535w^) - 5^^(859^^ 
-24765^2^ + 3333tu'^ - 2489m^) + 6st^{ll39t^ + 966t^u - 5034^1^^ - 17741^^) 
-24t^(t + u){8t^ + 355tu + 311u^)] + 8Q^[3364s^ + s^(26777t + 5430m) + 2s^{21659t'^ 
+9239tu - 4422^2) + 2/(14975t^ + 13824t2M - 13008tM2 - 6814?/) + 2s''^(17155t^ 
+14097t\ - 4659t\2 - 2969tM-'^ - 1116m^) + s2t(3214t^ - 12410f\ - 43626t'^u^ + 15952^^2 
+2941m^) + 6st^{t + M)(471f'^ - 1622t\ + 3626tM=^ + 3676^^) + 18f'^(t + uf{29t'^ 
-126tu + 367^2)] - 4g6[1872s^ + 8s^(2931t - 898m) + ^(88047^2 - 24599tM 
-23750^2) + 6s^(12809t^ - 3958t2M - 1260Um2 - 1538^3) + 2s^(43800t^ - 15493i^M 
-32520t2M2 - 9636tM^ + 5567m^) + 2s3(20619t^ + 7796t^M - 11308t^M2 - 2At\^ 
-3997t?/ + 2834?/) + 3sH{t + M)(2800t^ - 4748t^M + 10520t^u^ - 5494^^^ - 4779m^) 
+ 18st'^{t + m)2(98/ - 276t\ + 4:53tu'^ - 1036?/) - 36t^(t + uf{88t^ 
-137tu + 75^2)] - 2g^[2656s9 + 8s®(857t + 1698n) - 6s^(10701t2 - 9155tM - 1344^2) 
-s''(90673t^ - 58884t2?i - 29715t?/ + 25426m^) - 2s^(49901t^ - 71901t\, - 4:1961t\^ 
+28391tM^ + 16518m^) - 6s^(15176f' - 7421/m - 15765t'^?/ + 452t^u^ + 7555t?/ 
+1587m^) - 2s^{t + M)(11653t^ - 22240/m + 232tV - 10754tV + 11156tM^ - 476m^) 
+sH{t + M)2(262t^ - 2336t^M - 1326t'^u^ + 6178tu^ - 14453m^) 

+6st^{t + M)^(592f'^ - 1782t2M + 209 W - 1142m^) - 6t^{t + u)\232t^ - 232tu + 67^^)] 
+g=^s[2720s^ + 128s^(155t + 27m) - 45^(889^^ - 9125iM + 2274m2) - 45^(12176*=^ 
-e>2m\ + 3138tM2 + 3383m=^) - s^(66007i^ - 107787t=^M + 17517tV + 50095iM=^ 
-10752m^) - 4s^(16169t^ - 3039U^m - 13023iV + 22817i2M=^ + 3506iM^ - 6942m^) 
-2s^{t + M)(14858i^ - 31832i^M + 17273iV + 23882i2M=^ + 1273iM^ - 8164m^) 
-As^{t + uf{2mt^ - imt^u - 2955iV + 3776t2M^ - 1528iM^ - 753m^) - st{t + uf{208t'^ 
+5392i^M - 10164iV + 9544iM^ - 5213m^) + 12t\{t + uy{232t'^ - 232tu + 67u'^)] 
-s'^[96s^ + 8s^(519i - 22m) + 45^(1395^^ + 1447iM - 260m2) - 6s^(625i^ + 358t\ 
+6mu'^ + 92m^) - s^(12099i^ - 2111i=^M + 23601t\'^ + 2067tu^ - 608m^) - s^(9513t^ 
-24944i^M + 1486iV + lU72t^u^ - 12b7btu^ + 1312m^) - 6s^(i + M)(596i^ - 3874t^M 
+168UV + mAfu^ - 2541iM^ + 648m^) - 2s'^{t + uf{AAAt^ - I9m% 
+3031iV + 2304i2M^ - 3716iM^ + 1444m^) - su{t + M)^(608i^ - 1576i^M + l8mW 
-1321iM^ + 704m^) - 3tu^{t + uf{232t^ - 232tu + 67m2)], 
A = 256g^S(49s2 - llsi + 243^^) - 128g^^s[425s=^ + s^{m + 68m) - 7st{3At + 163m) 
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+9i^(135i - Uu)] + 192g^2s[565s^ + s^{753t + 317u) - s^{624:t'^ + 227tu - 2Alu^) 
-st(13Af + UAltu - 643x^2) - 3t^{7nt^ + imu + 2391^^)] - 32g^°s[3523s^ 
+5s^(2243i + 348ii) + 6s^(265i^ + hmu - Ahlv^) + 2s'^{QZm^ + 9495^^^ + 5907^1*^ 
-1123ii^) - st(25922i^ - m%7t^u - 19623iw^ - 4183ii^) + 9t^(285i^ + IMQt^u 
+1670iii^ + hhQu^)] + 16g^s[2564s^ + 19s^(1168i - 2?,hu) + s^(10973i^ - Zimtu 
-13434^2) + 2s3(32611t^ - MSt^u - 13905iw^ - AQ2hu^) + 5^(3319^^ + bl22At\ 
-24192i^w2 - 33026^1*2 - 1321m^) - si(2599i^ - 34853i^'u + 9564^2^2 + 32050^1^^ 
+12817ii^) + Qt^{t + M)(180t^ - 21lt^u - imtu^ - 499^3)] + 24QS[880/ 
-10s^(594t - 455m) - 5^(5508^^ - 2,imtu - 3085^2) - s^(20143t3 - 24746t\ 
-5389tM2 + 6032m^) - s^(22347t^ - 9890t\ - 272m\^ + 9660tM^ + 8167m^) 
-s^(9t^ + 3014t^M - 33328t^M^ - ?,^7Qt\^ + 11363tM^ + 2720m^) - st(t + ii)(2780t^ 
+ l^m^u - U91t\,^ + 4661tu^ + 3803?/) - 3t^{t + uf{384:t^ - mAfu 
+415tM^ + 194^3)] - 4gS[4544s^ + 4s^(56t + 2091?7,) - 12s^(481t2 - 792tM 
+419^2) - s^(18620t^ - 40275t2M + lb?,mtu^ + 11117^3) - /(51008/ - lU^mt^u 
+lUA<9t^u^ + 36149tM^ - 11325m^) - 3s^(12373t^ + 1095/m - W92Qt^u^ + 30592t\^ 
+Qmtu^ - 7111m^) + s^{t + ?7,)(2432t^ + 20995t^M + 7298^^^^ - MUQfu^ + bbl^tu^ 
+7739m^) - st{t + M)^(8896t'^ - IQhSQt'u + 10257tV - 2351tM^ - 3821m^) 
+9t^(t + M)3(240t^ - 2ASt^u + 68^2 + 25m^)] + 2Q'^.s^[m2s^ + 8s^(832t - 49m) 
+45^(1444^2 ^ 2143tM - 1509^2) - 2s^(4922t3 + 7689t^M - 2772tM2 + 2165m^) 
-5^(15173/ - 19913t^M + 60255t^M^ - 21787tM'^ - 920m^) - 4/(2242t^ - 13469t^n 
-hmit^u^ + l{}3hAt^u^ - 6857t?/ + 807?/) - 2s'^{t + M)(3892t^ - 2872/m + 775^^^^ 
+2434t2M3 _ 7345^^4 ^ 3106^^) + 2s{t + uf{lim^ + 3888/m - SlUt^u^ 
+2Smt\^ + mQtu^ - 1153m^) - tit + m)^(272/ + 2048t^M - 2A72t^u^ 

+9>mtu^ + 2mu% (A. 14) 
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(a) (b) 

Figure 1: Representative Feynman diagrams for the partonic subprocesses (a) i^e + — > 
e- + cc[n] + u and (b) Ue + s^e- + cc[n] + c, where n = ^S[^\^Pi^\^S^^\^S[^\^PP . In 
case (a), there is one more Feynman diagram, which is obtained by connecting the gluon 
hne with the -u-quark hne. In case (b), there are three more Feynman diagrams: two of 
type (a), where the d and u quarks are replaced by the s and c quarks, respectively, and 
one of type (b), where the gluon line is connected with the c-quark line to the right of the 
csW'^ vertex. 
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Figure 2: NRQCD (upper histogram) and CSM (lower histogram) predictions for the 
Pt distribution of prompt J/ip inclusive production in i/N CC DIS appropriate for the 
CHORUS experiment. 
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Figure 3: NRQCD (upper histogram) and CSM (lower histogram) predictions for the 
y distribution of prompt J/-^ inclusive production in vN CC DIS appropriate for the 
CHORUS experiment. 
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Figure 4: NRQCD (upper histogram) and CSM (lower histogram) predictions for the 
Pt distribution of prompt J/xjj inclusive production in ep CC DIS appropriate for the 
THERA experiment. 
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Figure 5: NRQCD (upper histogram) and CSM (lower histogram) predictions for the y 
distribution of prompt J/V' inclusive production in ep CC DIS appropriate for the THERA 
experiment. 
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